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Abstract. Local-type primordial non-Gaussianity couples statistics of the curvature perturbation 
( on vastly different physical scales. Because of this coupling, statistics (i.e. the polyspectra) of £ 
in our Hubble volume may not be representative of those in the larger universe - that is, they may 
be biased. The bias depends on the local background value of (, which includes contributions from 
all modes with wavelength k < Hq and is therefore enhanced if the entire post-inflationary patch is 
large compared with our Hubble volume. We study the bias to locally-measured statistics for general 
local-type non-Gaussianity. We consider three examples in detail: (i) the usual /nl, 9nl model, 
(ii) a strongly non-Gaussian model with £ ~ (q, and (iii) two-field non-Gaussian initial conditions. 
In each scenario one may generate statistics in a Hubble-size patch that are weakly Gaussian and 
consistent with observations despite the fact that the statistics in the larger, post-inflationary patch 
look very different. 
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1 Introduction 

A primary goal of observational cosmology is to learn about the physics of inflation. A key observable 
is the statistical distribution of the primordial curvature perturbation In single-field models of 
inflation the statistics of £ are inherited from quantum fluctuations in the inflaton itself [1-5]. If 
multiple light fields are present during inflation, such as in the curvaton scenario [6, 7] or modulated 
reheating [8, 9], fields other than the inflaton may generate £ and the relationship between the 
post-inflationary curvature and quantum fluctuations generated during inflation can be non-linear 
(see e.g. [10]). In these scenarios, the inflaton and curvaton field may obey Gaussian statistics while 
the observed curvature perturbation is non- linearly related to the Gaussian fluctuations, allowing 
for the phenomenological parameterization: 

C(x) = Cg(x) + 3 -f NL (£(x) - (C 2 G )) + ^9NL (CS(x) - 3«£>Cg(x)) + • . . (1.1) 

where Cg is a Gaussian random field and /nl, 9nl are constants specified by the particular in- 
flationary model [11-15]. The field C tn Eq. (1.1) obeys non-Gaussian statistics and the level of 
non-Gaussianity can be characterized by the products /tvla/ (Cg)' 9nl{Cg)- • ■ 
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Anisotropies in the cosmic microwave background (CMB), along with other cosmological datasets, 
provide stringent constraints on the statistics of £ within our Hubble volume. The variance of 
the primordial curvature perturbation is determined to be A^(fe) = A^ VMAP (k/k p i vo t) ns ~ 1 where 
&WMAP = 2 - 464 ± °- 072 x 10 ~ 9 > k pivot = 0.002/Mpc and n s = 0.9608 ± 0.0080 at 68% con- 
fidence [16]. The non-Gaussian parameters in Eq. (1.1) are bounded —10 < f^L < 74 and 
-12.34 x 10 5 < g NL < 15.58 x 10 5 at 95% confidence [17, 18] 1 . In this paper we ask the fol- 
lowing question: Suppose that inflation lasted sufficiently long that our Hubble volume is small 
compared to the entire volume generated during inflation. Are the statistics of the curvature per- 
turbation observed in our Hubble volume necessarily representative of the statistics of the curvature 
perturbation in the rest of the universe? 

The cosmological principle might make this question seem unnecessary, but (as has been noted 
by a number of authors [20-27]) even if the universe is statistically homogeneous and isotropic, 
mode-coupling introduced by non-linear terms like those in Eq. (1.1) correlates the statistics of ( 
on very different scales. The locally observed, smaller-than-Hubble-scale statistics depend on the 
unobservable, super-horizon modes of £g (i.e. long-wavelength modes of the Gaussian field that 
are nearly constant across our Hubble volume). Changes to the local statistics depend on a small 
parameter, 0(( sup er —horizon) *S 1> but can nevertheless be important. In particular, the background 
mode in our Hubble volume is roughly a sum over all modes with wavelength ir/k larger than c/Hq, 
so the variance of the total background fluctuations is larger than the locally observed variance in 
a single k mode by an amount dependent on the number of super-horizon e-folds. This situation, 
in which long-wavelength modes of ( bias the local statistics, is in stark contrast to the case where 
£ is Gaussian and different Fourier modes are strictly uncorrelated. For a Gaussian field, the local 
power spectrum may be randomly different from the globally averaged one, but is not systematically 
biased. 

The possibly biased nature of local statistics in non-Gaussian cosmologies due to long- wavelength 
modes has also been investigated recently by Nelson and Shandera [26] and Nurmi, Byrnes, and Tasi- 
nato [27] (which appeared while this paper was in preparation). Nelson and Shandera studied the 
scaling behavior and squeezed limits of n-point correlation functions computed by observers with 
access to a finite subvolume of the universe. They pointed out, for the first time, that some regions 
can be so biased that the local statistics appear only weakly non-Gaussian, even if the statistics 
in the larger universe are strongly non-Gaussian (a situation we consider in detail in §4). Nurmi, 
Byrnes, and Tasinato study possible relations between the local model parameters in Eq. (1.1) as 
measured by an observer in our Hubble volume. While these authors focus slightly more on the 
implications of a potential non-detection of primordial non-Gaussianity, there is significant overlap 
with our analysis in §3. 

In this paper we study the bias to the local curvature perturbation field ( for non-Gaussian 
initial conditions that can be written as a non-linear function of a Gaussian field £g that is local 
in configuration space. In §2 we introduce our notation and formulate the calculation of local 
statistics of Cng m terms of a short-long wavelength split of the Gaussian field Cg- I n §3 we study 

lr There are, of course, a number of theoretically-motivated forms of primordial non-Gaussianity that are not captured 
by Eq. (1.1) (see, for instance [19] and references therein). However, in this paper we restrict our attention to "local" 
non-Gaussianity which can couple Fourier modes of vastly different wavelengths. 
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the mapping between the global and local values of non-Gaussian parameters in the case when 
the statistics in the larger universe are weakly non-Gaussian as in Eq. (1.1) (see also [27]). In 
§4 we present full calculations of a somewhat counterintuitive example first discussed in [26] in 
which the local statistics of £ can appear to be nearly Gaussian while the global ones are strongly 
non-Gaussian. In §5 we study the relation between globally and locally determined statistics for 
two-field initial conditions. In §6 we summarize our results and discuss some possible implications 
for interpreting potential measurements of primordial non-Gaussianity. Appendix A contains a 
diagrammatic formulation of non-Gaussian statistics for a general local model, which is useful both 
computationally and conceptually for understanding the mapping between local and global statistics. 
For the brave at heart, calculations of global and local statistics for generic local- type non-Gaussian 
initial conditions are given in Appendices B and C and discussion of diagrammatic techniques in 
Fourier space is given in Appendix D. 



2 Local Statistics in a Subvolume 

We are interested in mode-coupling for the non- linear curvature perturbation £. Following [28], 
we let a be the locally defined scale factor and factor a into a spatially homogeneous piece and a 
perturbation, 

a = ae c(x) (2.1) 
where £ is a perturbation with volume average equal to zero over a volume Vl so that, 

a(t) = Iv l d ' x ln " and C(x) = In Q . (2.2) 

The locally observed perturbation with respect to some smaller volume Vs will be 



Cs(x) = In a - — / d 3 x In a (2.3) 
IVs 



Vs 



= lna-— / d 3 xlna-(— / d 3 xlna-— / d 3 xlna) (2.4) 
v l Jv l \Vs Jvs v l Jv L ) 

V v ' V v ' 

c(x) =aw 

So that, 

Cs(x) = C(x) - Cl(x) (2.5) 

i.e. a change in reference volume amounts to an additive shift in the perturbations. 

Suppose the non-Gaussian statistics in a volume Vl can be written as a local, non-linear 
transformation of a Gaussian field as in Eq. (1.1), or more generally, (ng = /(Cg) — (/(Cg))- The 
non-linearity of / in the variable Cg couples Fourier modes of (jvc ° n different scales. Given a non- 
Gaussian curvature field defined in a large volume Vl, how will the statistics appear to an observer 
with access to only a finite subvolume V5? 

From Eq. (2.3), the local value of ( in Vs is 

C(x)| s = C(x) - J ^j k -*W s (k)((k) (2.6) 
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Figure 1. Quantities in our Hubble volume are different from those in the larger universe by terms 
proportional to Cg.l ~ the long wavelength fluctuations that appear to be constant within our Hubble volume. 
Plotted is (Cg.l) 1 ^ 2 assuming that A G = 2.464 x 10~ 9 {k/k piv ) ns ~ 1 for two constant values of n s , along with 
two examples that include running. Note that (Cgl) 1 ^ 2 is larger than the 10 -5 amplitude in an individual 
Fourier mode because (£ G L ) is roughly a sum over (£ G (k)) for k < H . Throughout this paper we assume 
constant n s — 0.9608 [16] or ji s = 1, but keep in mind that the results for N > 20 are very sensitive to the 
(unknown) infrared behavior of A G (k) . 



where 



giving 



1 



'3„ „— tx k 



W s (k) = ^r I tFxe 
vs Jv s 



C(k)| s = C(k)(l-W s (fc)) 



(2.7) 



For simplicity we make the approximation that W$(k) is a top-hat function in Fourier space so that 
we can define 



Co. 



,5(x) = / 



(2tt) 3 



Cc(k)e 



ik-x 



and Cg,l(x) 



|k|<fc 



d 3 k 



Cc(k)e 



ikx 



(2.9) 



1/3 

where k* ~ 2tt/V s . Now the local Gaussian curvature can be split into short and long- wavelength 
pieces, 

Cg(x) = Cg,s(x) + Cg,l(x) (2.10) 
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with auto-correlations 



f d 3 k 

(C G) 5(x)C Gl s(x')) = J ^3^(fc)e !k - (x - x,) (2.11) 
f d 3 k 

l\k\<k, (27T) 3 



d 3 k 

l\k\>k 

(Cg,l(x)Cg,l(x')> = / £^Pa(k)e°*W (2.12) 



and vanishing cross-correlation 2 

(Cg,l(x)Cg,s(x'))=0. (2.13) 
We've defined the power spectrum of the Gaussian field as 

(C G (k)C G (kO) = (27rfP G (k)5 D (k + k'). (2.14) 

In what follows we assume, 

(&,s)s*(&,s) while <Cg,l>s«Cg,l (2-15) 

where ()g indicates averages over the volume Vs and () indicates an average over the entire volume Vl- 
That is, we assume that the locally measured, small-scale Gaussian power spectrum is representative 
of the globally defined one and that the variation in long-wavelength modes (k^ < k*) across the 
volume Vs is negligible. 

In this limit an observer in the small volume Vs is unable to distinguish between Cg,l and the 
background. However, the nonlinear coupling of Cg,l to short wavelength modes Cg,s will cause the 
local, small-scale non-Gaussian statistics of Cng to differ from the global ones in a way that depends 
on the local value of fa,L ~ that is, the local statistics are biased by (g,l- 

The typical size of the bias is characterized by the variance of long- wavelength fluctuations, 

, 2 . f k * dk 



<Cg,l> = J a Y A G(k) (2.16) 

where A G (k) = j^ysk 3 'Pc(k) and A ~ 2ir/V^ 3 , the infrared cutoff corresponding to the larger 
volume where the perturbations are set up. In the example calculations and plots we take Vs to be 
our Hubble volume ~ H^ 3 , but the expressions in the paper are completely general. The general 
expressions may be relevant for making comparisons between theory and particular observables 
measured in a volume smaller than our Hubble volume. 

Letting the power spectrum for fa be a power law, din A G / dink = n s — 1 with n s — 1 = const. 
gives closed-form expressions for the variance of long- wavelength modes 

<CG,L>| ns=1 = A G iV (2.17) 

(I _ e -(n 3 -l)N\ 

(Cg^U = A G W 1 - _ - ) (2.18) 



2 The more realistic assumption of a top-hat window function in real space with radius V^ 3 generates ~ Ag(fc ~ 
2.5/Vg //3 )/4 corrections to Eq. (2.13), where kVg^ 3 ~ 2.5 is the peak of Ws(k)(l — Ws(k)) for a top-hat in real space. 
For n a = 1 this gives (Cg,l(x)Cg,s(x'))/<Cg,l(x)Cg,i,(x')) ~ 1/(4JV), and (Cg,l(x)Cg,s(x'))/(Cg,s(x)Cg,s(x')) ~ 
-l/(41n(|x-x'|iT )) for |x-x'|/V s 1/3 < 10. 
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Dictionary of Frequently Used Symbols 



Quantity 


Definition 


Cng, PNc(k) 


The non-Gaussian curvature perturbation and its power 
spectrum 


Co, P G {k) 


A Gaussian random field used to generate £atg and its power 
spectrum 


A|(fc) 


The variance of fluctuations of the field X, A^(k) = 

k 3 P X (k)/(27T 2 ) 


Vr 


The volume over which Ovg is defined, e.g. the entire post- 
inflationary patch 


Vq, Wq 

* J 5 v ' J 


A subvolume of Vl and the corresponding window function, 
Wg(k) = Jy^ (i 3 xe _ix ' k . For most of this paper we take Vs 
to be our Hubble volume 


Cd <?, Cc J 


The short and long wavelength components of C,g- See 
Eq. (2.11), Eq. (2.12) 


N 


The number of super-horizon e- folds, |ln(V^/V5) 


A 2 


The amplitude of scalar-perturbations at k p i v = 0.002/ Mpc 
as measured by WMAP [16] 


n s 


The scalar spectral index n s — 1 = (fin A 2 /dink 


i * h « ~n n 

AC 5 U 


The spatial curvature and Hubble scale today 


TNT, Qm t « tnt + hivr 


Non-Gaussian parameters given in Eq. (3.1) (or as defined 
by the squeezed limits of the bispectrum and trispectrum in 
Eq. (B.19), Eq. (B.21) and Eq (B.20)) 


lie 


The value of the quantity X measured in Vs 




In §3 and §4 we rewrite GvgIs i n terms of Xg( x ) oc Cg( x ) 


0"G, </>G, £ 


Gaussian random fields we use to define (ng i n the two- 
field example in §5, and the ratio of their power spectra, 

= P^/Pa 



where iV = ln(ifo/A) is the number of super- horizon e-folds from the start of inflation to the 
time when the comoving scale of the observable universe crossed into the horizon. In Figure 1 we 
plot ±(^ L ) 1 / 2 , the typical amplitude of the unobservable background mode, as a function of the 
number of super- horizon e-folds N. The power spectrum of £g is of course unknown for k < Hq, 
but as a starting point we consider constant n s as in Eq. (2.17)-Eq. (2.18). As can be seen from 
Figure 1, even the modest red-tilt that is currently favored (n s = 0.9608) dramatically increases the 
typical amplitude of super- horizon fluctuations relative to that for a flat spectrum n s = 1. This 
difference becomes significant for iV > —{n s — 1) _1 ~ 25 e-folds - precisely when 0((n s — l) 2 ) 
contributions to the running are expected to change n s — 1 by order unity [29]. The specific shapes 
of (Cgl) 1/2 Pitted in Fi gure 1 should therefore be interpreted with caution, particularly for N 3> 
l/(n s — 1). For reference, we also plot examples of (Cgl) with running spectral indices given by 
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Long wavelength modes of (jvg wm & lso contribute to the mean spatial curvature measured 
within our Hubble volume 3 

—2 f r/ 3 k 

V k \ s = ^pj j^k 2 ( NG (k)W s (k) . (2.19) 

In a given subvolume V5, knowing the value of the background mode (g,l, is insufficient to specify 
fifcL. Nevertheless, we can estimate the typical amplitude of O^L in the scenarios we consider 



— , ^k^ c (k)\Ws(k)\ 2 = 4A ' W (l - e -(«.+»)^ « 4 f°™ (2.20) 
9i/ 4 7 fc fV ;I M ;I 9(n, + 3)V / 9(n s + 3) 1 ; 

where in the final ~ we have assumed that the power spectrum is not too red (e.g. for n s + 3 > 1). 
So, the dominant contributions to come from modes with A; ~ -ffo and, in contrast to Eq. (2.17)- 
Eq. (2.18), there is no enhancement from N S> 1. We therefore ignore constraints on C,G,L coming 
from constraints in because only the first few modes outside the horizon lead to spatial curvature 
and in fact we are working in the limit that Cg,l is independent of x within our Hubble volume 
anyway (e.g. Eq. (2.15)). However, it would be interesting to revisit these constraints and their 
implication for the bias of local statistics, particularly if local-type non-Gaussianity is detected. 



3 Example I: Weakly Non-Gaussian Initial Conditions 

In this section we imagine that the statistics in the larger volume Vl can be described by the usual 
local ansatz 4 

Gvg(x) = Cg(x) + \f NL (C G (x) - (C 2 G )) + ^g NL (Cg(x) - 3(C g )Cg(x)) (3.1) 
+ lk hNL (Cg-6(Cg)Cg(x) + 3(C g ) 2 ) 
where the non-zero coefficients satisfy 

1 » fNL^h) » 9 nl(Cg) » W(Cg» 3/2 • ( 3 - 2 ) 

The equation above is the definition of weak non-Gaussianity for this model. Single-source non- 
Gaussian models with coefficients with this scaling will generate non-Gaussian polyspectra that scale 
as (C]vg) c ~ (Ajfg)" -1 , where Ajvc is the observed variance (for further discussion see Appendix 
B). We can then apply the condition in Eq. (3.2) to require that the power spectrum of (jvg ° n 

3 Here, we are using the scalar curvature on spatial hypersurfaces 7?' 3 ' = — 4V 2 £(x), and taking ttk\ s = 
-Affo f H -3 d 3 x7? (3) (x), however see [30-35] for more detailed discussions of constraints on ((k) contributions to 
Qk as measured in our Hubble volume. 

4 Here we subtract the 3<7jvl(Cg)Cg(x) so that the power spectrum is unaffected at linear order in gNL and subtract 
6/ijv£(CgKg( x ) so that the bispectrum is unchanged at linear order in /ijvl- This helps to isolate how each coefficient 
changes the statistics of £, particularly in the case where a lower-order coefficient is vanishing (e.g. /jvl = 0, but 
hNL 7^ 0). See Appendix B for general expressions relating the coefficients in Eq. (3.1) and series coefficients in a 
general local map between Cg(x) and £]vg(x) 
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Figure 2. The observed amplitude of scalar fluctuations in our Hubble volume systematically differs from 
the average value in the entire universe by a fractional amount 12/5/jvlCg.l where Cg,l is the (unobservable) 
background mode in our Hubble patch. Plotted is an estimate of the probability distribution for the true value 
of A 2 , given the locally observed value Awmap = 2.464 x 10 -9 for two values of /nlIs an d two values of AT, 
the number of super-horizon e-folds. An observer in our Hubble volume cannot measure N, and therefore is 
unable to determine which probability distribution correctly describes our universe. The vertical dashed lines 
show the 68% confidence interval (±0.072 x 10~ 9 ) on A^ MAP from the eCMB+BAO+HO dataset [16]. 



CMB scales agrees with the power spectrum of to some accuracy, that is we could require that 
the 0(fx L ) terms are not important. Note that depending on the shape of the power spectrum, this 

requirement may be much stronger condition than requiring that f^L \J ^wMAp(k), 9nl^-\ymap(^)^ 

hNL^\YMAp(^) < 1 on CMB scales. We have checked that the examples plotted in Figures 2 and 
3 satisfy Eq. (3.2) for the assumed n s = const, power spectra. 

In the larger volume, Vl the field Cng( x ) given in Eq. (3.1) has power spectrum 

P NG (k) = P G (k){l + 0((&))) , (3.3) 

and the bispectrum and trispectrum are characterized by the coefficients /nl, 9nl and t^l = 
(|/atl) 2 up to corrections {{($)). 

An observer in a finite region V$ with background field value Cg,l wm see local statistics 
described by 

Cng\ s = XG(x) + I f NL \ s (xg(x) - (Xg)) + 4 9NL\s [Xg - 3(x 2 >XG(x)) + . . . (3.4) 
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Figure 3. The observed amplitude of Jnl m our Hubble volume, Jnl\s, differs from the average value in 
entire universe depending on C,g,l, Inl and gNL- Plotted is an estimate of the probability distribution for the 
true value of Jnl given A^ My4p = 2.464 x 10~ 9 and fixed gNL — 10 3 (left panel) or g NL — 10 5 (right panel). 
For Jnl =30 and gNL = 1000 there is a partial cancellation between the gNL and f NL terms in Eq. (3.8) 
causing the distribution to be narrower than in the fNL = 10, 50 cases. The width of the distributions depends 
on the RMS value of Cg,l- We assume n s — const. — 0.9608, which in comparison to n s = 1, gives a difference 
in (Cg,lV /2 of ~ 10% at TV = 1 and a factor of - 3.5 by N = 100. 



where we've defined 



V;(X!- (l+|/^CG,L + 0(C^ L ))«, ( ;.s (3.5) 



which we require to give the locally observed power spectrum 

< XG (k) XG (k')) = (2^) 3 fe(k + k')P x (k) where P x (k) = ^mapI ' . (3.6) 

The local power spectrum P x is related to the globally defined one Pjvg(^) through 

P NG (k)\ s = P x (k) = (l + ~~ JnlCg,l + 0(C g ,l)) Pc(k) 
and the locally observed non-Gaussian parameters are 



(3.7) 



9 12 

fNL\s = fNL + ^9Nl(g,L - InlCg,L + 0(Cg,l) ( 3 - 8 ) 

12 18 

9nl\s = 9nl + -^h-NiCcL — ~jr Inl9nlCg,l + 0(Cg,l) • 

Eq. (3.7)-Eq. (3.9) show that the connected n + l-point functions of Cng adjust the n-point functions 
of Cng Is by terms 0(Cg,l) and cause the locally observed statistics to differ from the global ones. 
For a strictly Jnl model (i.e. gNL, huL,- • • = 0) with /ml > 0, a positive background fluctuation 
Cg,l boosts the local power relative to the local bispectrum, the net effect is to make the local 
statistics appear more Gaussian than they are in the larger volume Vi (i.e. Jnl\s < Inl an d 
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Inl\s y < /iVLV^A 2 ). Negative /nl or background fluctuations will, of course, have the opposite 
effect. On the other hand if g^L ~ /jvx> then leakage from the trispectrum into the bispectrum 
can compensate and the local /nl value can be representative of the globally defined one. The 
cancellation between ^gNL and Tp/jvL in Eq. (3.8) is precisely what happens in the curvaton model 
when the curvaton dominates the energy density of the universe at the time of decay [36]. However, 
the level of non-Gaussianity as quantified by /jvlVA 2 is still adjusted. 

If Jnl 7^ 0, the measured value of the scalar power spectrum in our Hubble volume A\ymap 
differs from the average value in the larger universe Vl by an unknown amount ^/nlCg,l _ unknown 
because we don't know the values of Cg,l or /jvx- In Figure 2 we plot an estimate of the probability 
distribution for A 2 in Vl for fixed values of /jvx, assuming the observed value is Al^ MAP ^ and that 
Cg,l is drawn from a power law spectrum as in Eq. (2.18). Similarly, the local /jvx an d g^L values 
in Vl are related to the observed ones by amounts dependent on (g,l- In Figure 3 we plot estimates 
for the distribution of Jnl values in Vl assuming the locally observed power spectrum, and several 
possible values of /jvzJs> 9Nl\s- The probability distributions plotted in Figure 2 and Figure 3 are 
estimates of the probability distributions in that: (i) we don't allow all the observed parameters 
(Aq\ s , Inl\si Snl\s) to vary simultaneously and (ii) we neglect terms 0((q l ) m relating values 
of parameters /nl, 9nl measured in our Hubble volume to those in the larger universe Vl- A more 
realistic, but more involved calculation would be to calculate the posterior probability distribution 
of (A 2 ^, Jnl, 9nl) given the observed values A 2 ^, Jnl\si 9nl\$ along with their observational 
uncertainties, and the fact that Cg,l is Gaussian distributed. We are merely interested in illustrating 
the range of possibilities and leave a thorough exploration of parameters for another study. We 
further emphasize that observationally, we don't have observational access to JV - a parameter we 
have held fixed in Figure 2 and Figure 3. The variation in the probability distributions for different 
./V values should therefore be interpreted as an additional observational uncertainty. 

3.1 Implications for Model Builder 

When constructing a model of inflation, one typically specifies some set of fields relevant for inflation 
and the primordial fluctuations, as well as any interactions the fields may have. This guarantees 
the existence of an inflating solution and fluctuations and determines the possible shapes of the 
correlation functions. Adjustable parameters then allow the model builder to match the observed 
amplitude of fluctuations and to tune any non-Gaussianity to an amplitude consistent with obser- 
vational constraints. The length of slow-roll inflation may or may not be an independently tunable 
microphysical parameter. How should the model builder decide if a given set of microphysical pa- 
rameters gives rise to a significant number of Hubble volumes consistent with the one we see? In 
non-Gaussian models the necessarily statistical nature of making predictions from inflation for our 
Hubble volume becomes much more important, even for relatively short durations of inflation. 

To illustrate this point, consider a very simple (if unrealistic) model with only quadratic non- 
Gaussianity, /nl, and a constant spectral index. Expressing the amplitude of fluctuations and of 
non-Gaussianity in terms of the parameters of the large volume theory and the subsample bias gives 
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a sense of how the local statistics can differ from the global statistics: 



/nl|s = /nl (l - ^hL^fiCh) , (3.10) 
/nlA| 5 = /nlA - ||/nla/^> B 

where (C G ) is the correlation function at zero separation - a constant 5 . We've defined the bias, B, 
as in [26] so that it increases as N increases and for fixed N is larger for rarer fluctuations: 

In this example we restrict to weak non-Gaussianity, so by the condition in Eq. (3.2), /nlv/ (Cg) ^ 1- 
Since the long wavelength modes are only a fraction of the total modes contributing to (Cg*)' f° r 
a scale invariant spectrum the bias is also less than one in magnitude except for extremely rare 
fluctuations, \B < 1|. 

Notice that for non-Gaussian inflation models, matching parameters in the theory to agree 
exactly with our local observations makes sense only if the number of e-folds in the model is not too 
large. One way of visualizing this criteria is plotted in Figure 4. If the number of e-folds in the theory 
is larger, the parameters should not be matched identically to what we observe on CMB scales. In 
that case sub-volumes that have statistics identical to the parent will be rare, and so our observed 
universe will not be the typical outcome of those models. Finally, we note that because both the 
amplitude of fluctuations and the value of /jvx are changing in typical subsamples as we look on 
different scales, it is useful to plot the quantity that shows how non-Gaussian the subsamples are on 
average. The relative amplitude of non-Gaussianity in the subvolume to that in the large volume is 
shown in Figure 4. Note that for positive /nl, an overdensity Qg,l > causes the non-Gaussianity 
to be smaller in the small volume. Similarly, an underdensity Cg,l < causes the non-Gaussianity 
to be larger. 

4 Example II: Strongly Non-Gaussian Initial Conditions 

Suppose the non-Gaussian curvature perturbation, Cng, i n the larger volume Vl is given by 

Ov G (x) = C£(x) - (C P G ) (4.1) 

where p is a positive integer > 1. This field has statistics that are not accurately characterized by 
an expansion of the form Eq. (1.1), in particular the polyspectra have a different shape and scale 
dependence from the local shapes given in Eq. (B.19), Eq. (B.21), and Eq. (B.20). Nevertheless, in 
the squeezed limits that observationally define JnLi 9nl, and t^l one finds 

/iVi-^i^, STVL-^p, r^vL-^p for p even (4.2) 



5 One might worry that we are scaling quantities by a loop factor (Cg), which is dependent on the power spectrum 
over the entire range of scales (and, without a cutoff is formally divergent for a scale invariant power spectrum). 
However, (£ G ) is merely a placeholder and the actual value cancels when calculating observed quantities - our results 
do not depend on the unknown UV behavior of the power spectrum Pg- 
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Figure 4. Left: The degree of variation in non-Gaussianity in subvolumes as quantified by the variance 
Z%g = (( ifN f Ly JS- lS ) 2 ) - 1 = (lfNL) 2 (C 2 G )(B 2 ) = {lfNL) 2 (C 2 G , L )- The horizontal, dotted line indicates 
Y>ng — 0.15, roughly when uncertainty in f^L due to super-horizon correlations becomes comparable to the 
expected error on Jml from Planck (we assume ±5), if /nl = 40 in Vl- Right: The fractional change in the 
level of non-Gaussianity, vs - number of super-horizon e-folds N. The upper and lower curves 

correspond to £g,l = (Cg.l) 1 ^ 2 and £g,l = — (Cg^l) 1 ^ 2 , respectively. In these figures we've assumed gNL, 
h NL . . . = 0. 



Jnl = , 9nl ~ , r NL ~ "-TT— for p odd . (4.3) 

In contrast to the weakly non-Gaussian case in §3, this field has 1 ~ Inl\J ^-%g ~ 9nlA 2 ng ~ 

tnlA- 2 ng , where A^ G ~ (Ccf) is the observed variance. In general, the /nl, 9nl and tnl will also 
be scale-dependent functions of k s , ki, the long and short- wavelengths used to take the squeezed 
limits in Eq. (B.19), Eq. (B.21), and Eq. (B.20). For a more thorough discussion of weak and strong 
local non-Gaussianity, see Appendix B. 

Consider the local statistics of (jvc in a subvolume of size Vs- The local non-Gaussian curvature 
can be written in terms of short and long wavelength modes of £c as in Eq. (2.11), Eq. (2.12) as 

Gv G (x)| s = j*S£c G ,s(x) + ^^yyCcS (C G ,s(x) - <C G ,s» + ^^y& G ,s + ■■■ (4-4) 

where we have suggestively ordered the series with the term linear in (^g,s first. Now, if we happen 
to be considering a small volume Vs with a background fluctuation satisfying 

(G^^JiCeJ, (4.5) 
then to a local observer CngIs given in Eq. (4.4) appears to be a field described by a weakly non- 
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Figure 5. Here we suppose that the curvature perturbation in Vj, is given by Cng( x ) = Cg( x ) p — ((g) i an< ^ 
ask whether the curvature perturbations in our Hubble patch could appear to be described by the weakly 



non-Gaussian series (e.g. Eq. (4.6)). We require a large local background fluctuation £g,l >• J{Cg s) an< ^ 
consistency with the observed power spectrum A^ MAP . Plotted are the corresponding values of Jnl (left 
panels) and gNL (right panels) for different values of the power law index p and Cg,l- Cg,l = \J {Cg l) (solid 

lines), (g,l = 3y \C, G L ) (dashed lines) and (g.l = §J{Q},l) (dotted lines). The upper row uses n s = 0.9608 
for A G , the lower row uses n s = 1. The bend in the plots for n s = 0.9608 occurs at N(n s — l)/p ~ 1. 



Gaussian expansion of the form 6 

Gvg(x)| 5 = XG(x) + I f NL \ s (Xg(x) - (xl)) + ^ 9nl\ s (x G (x) - 3 X g(x)(xg)) + • • • (4-6) 

6 Here we're assuming p < 10, say, or small enough that the binomial coefficients k] ^p!_k)\ don't spoil the smallness 
of the quantity fc!(p P i fc) , JiCaJ/Cod when J((gJ/Cg,i < 1 as given in Eq. (4.4). 
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where 



Xg(x) = K g ^Cg,5(x) + [ ^> ] , (4.7) 



G,L 



(the ( ((q s) Kg l ) ^ s because we have subtracted 3 gNL\$ Xg(Xg) fr° m the linear term in Eq. (4.6)) 
and 



3 , , _ P- 1 ,,J^hA „„, 9 , _(p-l)(p-2) , „AC| 



G,S/ 




Now, < 1 so the field in Eq. (4.6) should have large local non-Gaussianity. However, it is 
possible for the local statistics to appear only weakly non-Gaussian, i.e. 

f*L\ S J^>9NL\ S {&),•• ( 4 - 9 ) 

on top of sufficiently large background fluctuations. Taking Vs to be our Hubble volume and as- 
suming that Cg has a power-law spectrum with constant spectral index as in Eq. (2.18) the criterion 
given in Eq. (4.5) for observing weak non-Gaussianity can be written 

1 (4.10) 

-(n.-X)N for Us = const - ^ 1 ( 4 - n ) 

where as, before N is the number of super-horizon e-folds and we have introduced N$, the number of 
sub-horizon e-folds. For reference, N$ ~ 60 g ivgs \j ~ — _ / — —\)n rN "' for TV rs -' 1 5 Enid \l ~ — _ / — —\)N r * J 
0.05 by N ~ 150 when n s = .9608. 

Equations (4.10), (4.11) show that even for strongly non-Gaussian statistics of Cng i n Vl, the 
statistics in Vs appear weakly Gaussian on top of very rare background fluctuations ((g,l/((g l) 1//2 }t 
1). But for very large N the statistics in Vs appear weakly Gaussian even for typical values of Cg,l- in 
regions where Eq. (4.5) is satisfied, the possible values of /nl\s m subvolumes depends qualitatively 
on the sign of p: for even p, /nl\s > in all subvolumes, whereas for odd p the sign of Cq l ls 
significant and Jnl\s can be negative. On the other hand, gNL\s > f° r all values of p. 

Now we ask what the restrictions on p, A G , and Cg,l are in order to generate a curvature 
perturbation as in Eq. (4.6) that satisfies the observational constraints on the power spectrum, f^L 
and gNL in our Hubble volume. If we fix the ratio Cg,l/ \J '(Cq l) ( w bich is a measure of the rarity of 
our Hubble patch), the index p, and the observed level of fluctuations A?. = A^ MAp , Eq. (4.7) and 
Eq. (4.8) allows us to solve for the variance of fluctuations in the (unobservable) background field Cg 
along with the observed values of Jnl arid g^L as a function of N. The results are plotted in Figure 
5. We see that current constraints on the observed level of non-Gaussianity are indeed compatible 
with a scenario in which our Hubble patch is a biased subsample of a larger universe with strongly 
non-Gaussian initial curvature perturbations Ovg( x ) = Cg( x ) ~~ (Cg)- 
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5 Example III: Two-field Initial Conditions 

In this section we consider initial conditions inspired by a version of the curvaton model [6, 7, 36, 37] 
in which perturbations from both the inflaton <j> and the curvaton a are responsible for generating 
C (see e.g. [38-41]). In this "inflaton-curvaton" scenario, the curvature perturbation in the larger 
volume is given by 

Gvg(x) = <fc(x) + <r G (x) + " <<4» (5-1) 

We assume that 4>g(x) and o"g( x ) are statistically independent (i.e. (0(k)cr(k')) = 0), Gaussian 
random fields with proportional power spectra 

C 2 ^^- so that P NG (k) = P a (k)^l + e + ^f 2 N Mk^j (5.2) 

where I a (k) ~ is defined in Eq. (B.12) and for simplicity we assume ^ is a constant (however, 
see e.g. [40, 42, 43]). 

5.1 Case 1: a Is Weakly Non- Gaussian 

First, we make the usual assumption that the curvaton contributions to the curvature perturba- 
tion are only weakly non-Gaussian. That is, we assume that Jnl\/ <C 1. The non-Gaussian 
parameters that characterize the bispectrum and trispectrum of £ are 

Jnl (lf NL y 

f"L = (l + £2)2 ' TNL = (l+£2)3 ( 5 ' 3 ) 

and gNL = 0. In a subvolume, Vs a local observer will see statistics described by 

Png\ s = Png (l + (5.4) 

JNLls = (1 + ^2)2 I 1 + y £r^\JNLOG,L I (5.5) 
SWlIs = ( 5 - 6 ) 
1/ 



r ^ls = (i + ^)3 y- + y ^t^ ct g,l J • (5.7) 



In contrast to the case in §3, the local statistics are now modulated by long-wavelength modes of gq 
only, as opposed to fluctuations in the total curvature fluctuation £g,l = 4>G,L + o~g,l- To compare 
with §3, we rewrite Eq. (5.4) - Eq. (5.7) in terms of /nl, in Vl, 

Png\ s = Png U + ~ Jnl (l + £ 2 ) <tg,l\ (5.8) 
fNL\ s = f NL { 1 + f ( (i)2rjV / ;" 2/ ^ '| (1 + e 2 ) <^ (5.9) 



7VL 
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Now, the amount by which the power spectrum and Jnl vary from place to place is the same as in 
Eq. (3.7) and Eq. (3.8) with Cg,l ->■ (1 + £ 2 Vg,l- The typical size of the modulation of statistics in 
Vs is yfijl + £, 2 ) 2(J g l) = yj 0- + £ 2 )(Cg l)- S° for fixed /jvx, the typical modulation in the power 
spectrum is larger relative to the case where a single field, Cg> generates density perturbations. That 
is, a one-sigma fluctuation in uq,l generates a larger change than a one-sigma fluctuation in Cg,l, 
larger by a factor yl + £ 2 . The a field itself must be more non-Gaussian to maintain a fixed /jvl 
in the curvature as the power from <r decreases (£ increases). In Figure 6 we plot estimates for 
the distributions of A 2 and /jvx hi the total volume, given several values of £ 2 . As in §3 these are 
estimates of the probability distributions in that we (i) neglect terms O(Cgl) m relating values of 
parameters measured in our Hubble volume to those in the larger universe Vl and (ii) we don't allow 
the observed parameters (A^L, Jnl\s) to all vary simultaneously. We have again fixed the number 
of e-folds for illustrative purposes even though this is also an unobservable quantity. 

Note, that while the observed amplitude of the three-point and four-point functions are char- 
acterized by /tvl y^A 2 , ttvlA 2 , when one considers the entire series of correlation functions neither 
product alone quantifies the level of non-Gaussianity in the field (ng- The single quantity that 



controls the level of non-Gaussianity is 



i+^r (rather than f N L^c - (I+frp. 
tity is small, the series of cumulants is ordered and the amplitude of each consecutive cumulant is 
smaller by this factor. Each cumulant also has an extra factor of 1/(1 + £ 2 ) which does not affect 
their relative importance. In terms of the observed non-Gaussian parameters given in Eq. (5.3), this 
criterion for weak non-Gaussianity, ^_^^ c <C 1, is equivalent to requiring the kurtosis to be much 



When this quan- 



smaller than the skewness: tjv^A 2 - <C / 



TNL < 



V 



As in the single field case, we can ask how the total 



amplitude of non-Gaussianity differs in biased subvolumes: 



Jnl< 



where the bias here is defined as 



Inl\ 



QfNL(C G 



-B 



B 



(5.10) 



(5.11) 



The relationship between the amplitude of non-Gaussianity in Vs and Vl has the same structure as 
in the single field case, but the bias will generally be smaller (assuming the same total amplitude of 
fluctuations, y^A 2 .) since the fluctuating field contributes only part of the total power. Notice that 
when £ = this reduces to the single field expression, Eq.(3.11). 

5.2 Case 2: a Is Strongly Non-Gaussian 

Now we assume the perturbations coming from the curvaton are strongly non-Gaussian /nl^A 2 ~ 
1, but a subdominant contribution to the total curvature (£ 2 3> /^ L A 2 ~ 1). To understand 
the dependencies, it's helpful to define the 0(ffj L ) fractional change in the globally defined power 
spectrum in Vl 



e(k) = 2 



IInl) Uk) 



(5.12) 
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Figure 6. In the two-field "weak non-Gaussianity" case the amplitude of scalar perturbations, A 2 ^ and of 
/nl in our Hubble volume, /nl\si differ from the average value due to the background value of <Jg,l, while 
the total curvature is set by (Jg,l + 4>g,l- For fixed /nl, changes to the local statistics are typically larger by 
a factor of y 1 + £ 2 relative to the case in §3. Plotted are estimates of the probability distributions for A 2 
(left panel) and Jnl (right panel) in Vl, given the observed values in Vs for different values of £ 2 = P^/P a . 
Note that for £ 2 = 1, f NL \ s = f NL + 0(a 2 G L ). 



By assumption e ~ l/£ 2 and again I a ~ A 2 as defined in Eq. (B.12). In the larger box Vl the power 
spectrum is 

PN G (k) = ePa(k) (1 + e(fc)) . (5.13) 

Taking the squeezed and squashed limits in Eq. (B.19), Eq. (B.21), Eq. (B.20) gives scale-dependent 
non-Gaussian parameters: 

fNL(k) = -^-e(ki), g NL = 0, T NL {h) = (-f NL (h)\ -z (5.14) 

where k[ is the magnitude of the long- wavelength mode used to calculate the squeezed and squashed 
limits. The scale dependence of /nl, t nl is given by the function I a {k) in Eq. (B.12). For > A, 
where A is the infrared cutoff in A 2 (/c), the scale dependence of I a {k) is generally weak: for n s = 1, 
I a (k) ~2A 2 ln(A;/A). 

In this example, the field C is weakly Gaussian with hierarchical cumulants in that 1 S> 
Jnl\J ((g) ^ t nl(Cg)- However, in contrast to the weakly non-Gaussian, single-source case in §3 
where the cumulants scale as (Cng 1 )c/ (Cng) c ~ (£ng)> *he hierarchy of cumulants in this example 
scales as 

(/-n+l\ J Z/-2 U/2 l C 2 \l/2 

We have assumed that £ 2 >> 1, but depending on the relative magnitudes of £ and (Catg)' the 
higher-order cumulants may be more important relative to the lower order ones than in the examples 
considered in §5.1 and S3. 
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In a subvolume Vg, an observer will see a local power spectrum 



P N a(k)\ s = ePa(k) I l + e s (k) I l + 1 1 ( 5 - 16 ) 



and non-Gaussian parameters 

/jvl- . 2cr G,L \ ^ I p. ^ (f/jvrls) 



7-mI.v ■ — f * I 1 + ,. , ,„ , J ■ ;/.\ l\s ■ ■ • ~.\ /. ■ ■ — f — — - - ) 



where I as is Eq. (B.12) with P(k) -> P„{k)\l - W s (k)\ 2 and e s is Eq. (5.12) with -> J CTs . In this 
case the difference between the local and global statistics is more complicated: since I a ^ I (Ts the 
local statistics in Vs differ from those in Vl even if <jg,l = 0. Averaging over a G l will recover the 
parameters in Vl- 7 To see this, note that P a I a ~ Ia s Pa s + 2-P CTs (<7l 2 ) + I aL P aL so that I a {k) — > 
I as {k) + 2(<t l 2 > for scales k > Vg 1/3 , so that (e s (l + 2a 2 G ^ L /I <Ts )) -> e. 

To study the statistics in subvolumes V5, consider the following example: fix (|/tvl) 2 A 2 = 1, 
then we find expressions for the locally measured non-Gaussian parameters in terms of £ and the 
observed power spectrum A 2 ^ MAp , 

a V WMAP 

For simplicity we'll assume that n s = 1 so that we can use the analytic expression for I a (k) given 
in Eq. (B.13). Eq. (5.18) allows us to choose values of the ratio of inflaton to curvaton power that 
give particular values of Jnl, t~nl- We now choose £ such that fNL(kpiv) = 10 when averaged over 
the entire volume Vl- However, since the local value of Jnl depends on a G l , observers in a finite 
volume can easily measure fNL(k p i v ) 7^ 10. In Figure 7, we plot the probability distribution of 
possible observed values of /nl\s an d t nl\s- We see that unlike the cases considered in §3, §4, 
and §5.2 the probability distributions are extremely skewed and there is a large offset between the 
median and modes of the distribution of values of /nl\s an d t nl\s- 



6 Conclusions 

Local type non-Gaussianity couples the small-scale statistics measured by an observer restricted 
to a small volume Vs to the unobservable, long wavelength modes (g,l that are nearly constant 
across Vs- In this paper we have systematically calculated the relationship between local and global 

7 The average of the small-volume polyspectra over the long wavelength modes must recover the large-volume 
polyspectra. However, since the parameters /nl, 9nl etc. are ratios of quantities dependent on the random variable 
Cg,l, when terms non-linear in C,g,l are important the relationship between }nl\ s and /nl, say, is generally more 
complicated. In this example the non-Gaussian polyspectra are dependent on a, but power spectrum is dominated by 
the Gaussian field <j>, and therefore averaging over og.l doesn't change the denominator in the ratios used to define 
the non-Gaussian parameters, and the expressions for /jvl and tml are easily recovered from /nl\ s and tjvl| s by 
averaging over 0%^. 
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Figure 7. In the two-field "strong non-Gaussianity" case <j 2 g L modulates the local statistics leading to skewed 
probability distributions for the observed non-Gaussian parameters. The solid curves are the distributions 
for the observed values of Jnl\s (left panel) and tnl\ s (right panel) for different values of N - the number 
of super-horizon e-folds (in this plot we set n s = 1). In each case we've fixed (|/atl) 2 A 2 = 1 and chosen 
the values of Jnl and £ (see Eq. (5.1) - Eq. (5.2)) to produce /jvl(^ = 0.002/Mpc) = 10 when averaged 
throughout Vl (however, the corresponding t nl values are different). The values of f NL and t nl when 
averaged over Vl are indicated by the vertical dotted lines. 



statistical quantities (the power spectrum, bispectrum, and trispectrum) in models with local-type 
primordial non-Gaussianity. We demonstrate through explicit calculation, that broad classes of 
statistical distributions for the curvature perturbation in the larger universe Vl are consistent with 
nearly Gaussian statistics observed in our Hubble volume. This many-to-one nature of the mapping 
between statistics in Vl and Vs is potentially a challenge for using statistics measured in our Hubble 
volume to infer the statistics in the entire universe. The framework outlined in §2 and Appendix 
B is general, but we study three examples in detail: the usual local ansatz in §3, an example with 
strongly non-Gaussian initial conditions coming from a single field in §4, and finally a two-field 
example in §5. 

For the weakly non-Gaussian statistics for (jvg m §3 we find, in agreement with [27], that 
values of the non-Gaussian parameters /nl, 9nl consistent with current constraints can cause the 
statistics measured within our Hubble volume to differ from those in the larger universe, even for 
a modest number of super horizon e-folds (N ~ 0(10), say). This is illustrated in Figures 2 and 
3. These results are dependent on the unknown behavior of the curvature power spectrum on 
super-horizon scales (examples for different IR extrapolations of the power spectrum are plotted in 
Figure 1). Figures 2 and 3 assume that the power law spectrum on super-horizon scales remains 
unchanged out to k/Hg ~ e N , which may be false - the true behavior of A 2 (k) for k < Hq could 
increase or decrease the typical amplitude of Cg,l- Typical changes to the level of non-Gaussianity in 
Hubble-size subvolumes are plotted in Figure 4 for both n s = 0.9608 and n s = 1. In a universe with 
local non-Gaussianity, constraints on global statistics (and therefore inflationary parameters) from 
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observations in our Hubble patch are necessarily probabilistic because the locally observed power 
spectrum, bispectrum, trispectrum are dependent upon the unknown value of the random variable 
Cg,l- While this has been known for a long time in the context of slow-roll inflation [44-47], we 
have shown that for inflationary models with local non-Gaussianity - either strong non-Gaussianity 
or merely observable levels of non-Gaussianity - the probabilistic relationship between theory and 
observations is important. 

In §4 we consider strongly non-Gaussian statistics, Ctvg( x ) ~ Cg m the larger universe Vl that 
can appear only weakly non-Gaussian on sufficiently large background fluctuations Cg,l [26]. We 
determine the restrictions on p, A^ g and Cg,l to produce statistics consistent with observations in 
our Hubble volume. The main results are illustrated in Figure 5. We see that for sufficiently large 
N, typical subvolumes (e.g. (g,l/ \J '(Cg i) ~ ^' corresponding to ~ 30% of Hubble-sized patches 
in the universe) will have statistics consistent with constraints on parameters in the weakly non- 
Gaussian ansatz Eq. (1.1), even if the curvature perturbation in the rest of the universe is strongly 
non-Gaussian. At very large N, this can be true of the vast majority of subsamples (not only la 
and higher fluctuations) depending on the infrared behavior of the power spectrum. In this sense 
weakly non-Gaussian statistics may be considered 'natural' as discussed in [26]. While §4 focuses on 
initial conditions that are a single power law C p the qualitative results should hold for more general 
forms of strongly non-Gaussian initial conditions and we provide a framework for these calculations 
in Appendices A and B. 

Finally, in §5 we consider an example in which the initial curvature perturbation is given by 
a sum of two uncorrelated fields, one Gaussian <p and one non-Gaussian with a quadratic coupling, 
o = o~g( x ) + /atz,(0g( x ) ~~ (°g))- Initial conditions of this type are consistent with observations for a 
range of values of /jvl and/or £ 2 = P^jP a . The qualitative difference between this scenario and those 
in §3 and §4 is that the field that modulates the local statistics, o~G t L, is only partially correlated with 
the total curvature perturbation C = 4> + o~. There is therefore greater freedom in finding statistics 
in Vl that map to weakly Gaussian statistics in Vg- For a only weakly non-Gaussian, the results 
are similar to those in §3 but, for fixed Jnl arid A 2 , the typical size of the modulation in local 
statistics is enhanced by a factor yj 1 + P$ / P a (see Figure 6). On the other hand, if f NL ((r G ) 1/2 ~ 1, 
the results are qualitatively different: the probability distribution for observed values of Jnl, tnl 
is highly skewed (see Figure 7). 

The calculations in this paper are completely statistical: we do not attempt to give a dynamical 
model that generates the examples of statistics in Vl we have considered, nor do we attempt to 
understand how differences between local and global statistics alter inferences about particular 
inflationary scenarios. The parametric forms of initial conditions we have considered in §3, §4, and 
§5 are simplified examples of initial conditions that can arise in the curvaton, or modulated reheating 
scenarios, but we have assumed that the parameters (A^ and the coefficients of the non-linear terms, 
for instance) can be freely adjusted to tune the statistics in Vs, which is not necessarily the case. 

Throughout this paper we have made the simplifying assumption that the background mode 
Cg,l is precisely constant across our Hubble volume (see Eq. (2.15)). In reality slight variations in 
Cg,l from modes with wavelengths not too much larger than c/Hq, and correlations between these 
variations and the small-scale statistics of C may be detectable [30-35, 48]. While the assumptions 
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we have made in Eq. (2.15) should be sufficiently precise for sub-horizon scales k 3> Hq, and 
contributions to Cg L horn k <C Hq, it would be interesting to explore the potentially observable 
corrections for k ~ Hq. 
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A Diagrammatic Representations of n-point Functions 

We want to calculate ra-point correlation functions of the non-Gaussian field Ctvg defined by 

Gvg(x) = /(Cg(x)) - (/(Cg)) • (A.l) 

The n-point functions of Cng can be written entirely in terms of two-point functions of Cg an d the 
derivatives of /. However, the expressions quickly get messy so it's helpful to use connected diagrams 
to keep track of the terms (see also [25, 49]). 

Dictionary of Diagrams: 

In this paper a line segment connecting two points 1 and 2 represents the real-space correlation 
function between the Gaussian fields Cg at two spatial points xi and X2. 

<Cg(xi)Cg(x 2 )> = 19 92 (A.2) 

while a double line segment indicates the square of the Gaussian correlation function 

(Cg(xi)Cg(x 2 )) 2 = 1 ' 2 (A.3) 

and vertices with multiple line segments indicate products of correlation functions connected to 
different points 3 



(Cg(xi)Cg(x 2 ))(Cg(xi)Cg(x 3 )) = 




(A.4) 



2 
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Circles represent Cg( x ) contracted with itself which is independent of x. For instance, we can write 



(Cg(x 1 )Cg(x 2 ))(Cg(xi) 2 ) 




2 



(A.5) 



O 



so it doesn't matter which vertex a loop is connected to. 



The Two-point Function: 

The two point correlation function of the non-Gaussian field Ovg is given by 



(Catg(xi)6vg(x2)) 



(/ (D)2 + f Wf 



(3) 



o 



+ 



V 



(A.6) 



+ 



2 (\ 



+...) 



+ 



where the ... in each parenthesis represent terms with higher-order loop contributions, and the 
in the last line indicate terms with additional shapes (but still 0(C G ))- 



The Three-point Function: 

The three point function of (ng contains five types of terms at 0(Cg) : At this order the terms can 
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be grouped into three distinct three-point function shapes: 

((ng(xi)(ng{x2)(ng(x3))c = 




2 



2 2 
3 

+ ((/ (2) ) 3 + -..) 1 

2 

+ ... 

The Connected Four-point Function: 

Finally, we compute the connected four-point function up to 0(Cg) : 
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( Cn G (xi ) (nG (x 2 ) CtV G (X3 ) CNG (X4 ) ) < 



(/ (i) ) 3 / (3) + f i?(/ (i) ) 2 (/ (3) ) 2 + J(/ (i) ) 3 / (5) lO 



V 2 



+ 3 perm. 



(A.8) 



+ 



\ 

(/ (l) )2(/ (2) )2 + ( / (l)( / (2))2 / (3) + ( / (D)2 / ( 2 ) / (4)) 

o 



V 2 3 2 



n-x 



+ 5 perm. 



4 1 



+ (/ (2) r 



.1 



+ 2 perm 
4 / 



+ ±/ (1) (/ (2) ) 2 / (3) 1^ 



3 2 



4 1 



+ 11 perm. 



+ ^(/ (1) ) 2 (/ (3) ) 2 



3 2 
+ 

4 1 




+ 5 perm. J 



+ / (D (/ (2) ) 2 / (3) 



+ 11 perm. J 



+ ^(/ (1) ) 2 / (2 V (4) ( 




+ 11 perm. 
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B Statistics in V L 

Scaling of the Non-Gaussian Cumulants: 

Suppose that the statistics of the curvature perturbation can be written as a non-linear transforma- 
tion of a Gaussian field that is local in real space: 

Cng = /(Cg(x)) - </(Cg(x)> . (B.l) 
We can calculate the statistics of Cng m terms of (Cq) <C 1 and derivatives of / 



(B.2) 

Cg=o 



For the moment, we'll ignore the shape dependence of the n-point functions of Cng and just consider 
the scaling of the cumulants in terms of (Cg) (we consider the shape dependence in §B). In terms of 
/ and Cg we have, 

(Cng) = (Cg) {(/ (1) ) 2 + \ ((/ (2) ) 2 + 2/ (1) / {3) ) (Cg) (B.3) 



+ A( 5(/(3))2 + 6/(2)/(4) + 3/(1)/(5) ) {(2g)2 + 



(Cng) = (Cl? {3(/ (1) ) 2 / (2) + ((/ (2) ) 3 + Qf {1) f {2) f {3) + ^(/ (1) ) 2 / (4) ) (Cg) + ■ ■ ■ } (B.4) 

(C^ G )c = (C G > 3 (4(/«) 3 / (3) + 12(/«) 2 (/( 2 )) 2 + (3(/( 2 )) 4 + 36/( 1 )(/( 2 )) 2 / {3) (B.5) 

+12(/ (l) ) 2 (/ (3) ) 2 + 18(/ (l) ) 2 / ( 2)/ (4) + 2(/ (l))3 / (5)) <£)+...) 

(C%g)c = (Cg) 4 (5(/ (1) )V (4) + 60(/«) 2 (/( 2 )) 3 + 60(/«) 3 / (2) / (3) + . . . ) (B.6) 



(Cng)c = &- 1 (n(/( 1 ))"- 1 / (n " 1) + n(n - l)(n - 2 )(/( 1 )f- 2 / (2) / (n ~ 2) 



+ 



where the subscript c indicates the connected part and . . . indicate terms higher-order in (Cg)- ^ n 
Eq. (B.3)-Eq.(B.7) we've kept a number of subleading (in (Cg)) terms in order to help illustrate the 
following points: 

• For / V and / (n) // (1) (C G ) (n-1)/2 « 1, cumulants scale as (Q G ) C ~ (Cg)™ _1 ~ (Clrc)"" 1 - 
So, higher cumulants are suppressed by powers of the observed variance. We refer to this type 
of statistics as weakly non-Gaussian. 

• If /' = the cumulants still scale with increasing powers of (Cg)> but the observed variance is 
(Cng) ~ (Cg) P w here p is the order of the first nonzero derivative of /, so the relative scaling 
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of each {Cng)° ^ s different 8 . In particular the suppression of higher cumulants can be weaker 
than in the /' 7^ case. 9 We refer to this type of statistics as strongly non-Gaussian. 

Further note that for series coefficients f {n+1) y/i^) <^ f {n) when /( n ) is non-zero, non-Gaussianity 
will first be evident in either the skewness or kurtosis. That is, (Cng) — (Ctvg) c and/or (C/vg) c > 
{Qg)c for n > 4. 

Shape and Scale-dependence of the n-point Functions: 

We now consider the shape and scale-dependence of the power spectrum, bispectrum, and trispec- 
trum of (ng m terms of the Gaussian field (a- The statistics of Cg are completely specified by the 
two point function: 

(( G (k)( G (k')) = (2ir) 3 5 D (k + k')P G (k) , A 2 G (k) = k 3 P G (k)/(27r 2 ) (B.7) 

where Su is the Dirac delta function. The power spectrum, bispectrum, and trispectrum of (jvg are 
defined through 

(GvG(ki)Ov G (k 2 )> = (2vr) 3 <5 D (ki + k 2 )iV G (fci) (B.8) 
(0v G (ki)GvG(k 2 )0vG(k 3 )> = (2vr) 3 <5 1) (k 1 + k 2 + k 3 )S(k 1 , k 2) k 3 ) (B.9) 
(C7VG(ki)Cjv G (k 2 )CiVG(k3)C7VG(k4)) c = (2Trf5 D (k 1 + k 2 + k 3 + k 4 )T(ki, k 2 , k 3 , k 4 ) (B.10) 

Fourier transforming the real-space diagrammatic expressions in Appendix A we find that the non- 
Gaussian power spectrum is given by 

PNG(k) = Pa(k) { ((/«) 2 + fMfW((h) +.■■)+ Ic G (k) Q(/ (2) ) 2 + ...)+...}, (B.ll) 

where is the two-point function at zero separation (a constant) and we have defined, 

1 f rfik' 

For a scale-invariant spectrum, n s = 1, this becomes 

^W = -A G {ln^^ + ^} and <&> = A&ln (^) (B.13) 

where k max is the UV cutoff of the power spectrum and A ~ 2-ir/L is the IR cutoff. For k < \Jek max h, 
(( G ) > I^ G (k), but for k > \/ek max A, (Cq) < ^Cg(^)- ^Cg(^) ^ s plotted in Figure 8 for n s = 1. 

8 What we really mean by f^ n ' — is ~ ((Cg, s )), our expansion parameter, so that the next order terms are 
comparable. 

9 To be more precise, the scaling of cumulants with the non-Gaussian variance depends on which terms are non- 
vanishing in fix) — Y^, n -^f x n . If the only non-vanishing term has an odd power n, then the cumulants scale with 
powers of {Cng), but generically higher terms scale with ~ (( G ) > (Catg) arl( l are l ess suppressed than in the /' 7^ 
case. 
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The bispectrum is given by 
S( kl ,k 2 ,k 3 ) = (P G (k 1 )P G (k 2 ) + 2 P erm.) ((/«) 2 / (2) + (/ {1) / (2) / (3) + ^(/ {1) ) 2 / (4) ) <&> + -) 

+ 1 ( J(ki, k3)PG(fci)P G (A;3) + 2 perm.)) ((/ (2) ) 3 + . . . ) 

/ f (1) f (2) f (3) \ 
+ {I{k 2 )P G {k l )P G {k 2 ) + 5 perm.) I J ^ J +...)+... 

where J(ki, k 2 ) ~ 0(A G ) is a function that depends on both the magnitudes of ki and k 2 , and the 
angle between them 

J(k '- k '» = p G ( fa) p G (b) / (^F Pc(|k ' - k,|)Pc(|k2 + k,|)Pc(t ' ) • (R14) 

Notice that in the squeezed limit J(k;,k s ) — > 1(h) + (D(kf/kg). In another squeezed limit: |k s | = 
\k' s \ = |k s + k' s |, J(k s , -k' s - k,) + J(k' s , -k s - k,) + J(k s ,k s + k' s ) -> 3/3 Cc (fc s ) where, 

/3 fG (fc s ) = ^A^fc,) y rf 3 x/- 4 (l + x 2 - 2x • k s )^" 2 (l + x 2 + 2x • k' s )^" 2 (B.15) 

So, in the squeezed limit the angular dependence vanishes and the J functions are just dependent 
on the magnitudes k s , h. 

Finally, the trispectrum is given by 

T(k!,k 2 ,k3,k 4 ) = (P G (fc 1 )P G (A; 2 )P G (fc3) + 3perm.) (B.16) 

(/ (1) )¥ 3) + (}(/ {1) ) 2 (/ (3) ) 2 + ^(/ {1) )¥ 5) ) (c 2 G ) + ■■■ 

+ (P G (h)P G (k 2 ) (P G (|ki + k 3 |) + P G (|ki + k 4 |)) + 5 perm.) 

((/W)2 ( /W)2 + ( / (D (/ (2))2 / ( 3 ) + (/W) 2 / (2) / (4) ) <C G > + • • 

+ (K(k 1; k 2 , k 3 )P G (fci)P G (|k 2 + k 3 |)P G (fc 3 ) + 2 perm.) ((/( 2 )) 4 + 
+ ((/ Co (fei) + / CG (fc 2 ))P G (fci)P G (fc 2 )(P G (|ki + k 3 |) + P G (|ki + k 4 |)) + 5perm.) 
l /( i) (/ ( 2 ) )2/ (3) + _ 

+ (P G (fc 1 )P G (A; 2 )P G (A;4)J(k 2 ,k4) + 11 perm.) (/ (1) (/ (2) ) 2 / (3) + ... 
+ (P G (h)P G (k 2 ) (J fG (|ki + k3|)P G (|ki + k 3 |) + / fG (|k 2 + k 3 |)P G (|k 2 + k 3 |)) 

+5 perm.) Q(/«) 2 (/( 3 )) 2 + . . . ) 

+ (I CG (k l )P G (k l )P G (k 2 )P G (k z ) + ll V evm.) (±(fW)*fWfW + ... 

where K ~ 0(A G ) and depends on the magnitudes of ki,k 2 ,k 3 and the relative angles between 
them 

r^p G (fc')p G (|ki + k'i))p G (|k 1 + k 2 + k'i)p G (|k 3 -k'i) 





Figure 8. Plotted is I^ G (k), the higher-order correction to polyspectra of Ovg given in Eq. (B.12). Shown 
here for a scale invariant spectrum Aq (e.g. Eq. (B.13)). 



In the squeezed limit needed to calculate g n l we find 

K(k s XM ^ (B.18) 

where in taking the limit we have fixed |k s | = |k' s | = |k s + k^|. 

It is perhaps more useful to consider the usual local parameters (see, e.g. [50]) : 

3, _ 1 g(kj,k a ,-k f -k 8 ) 

5 /M = 4ft PffiPffi (R19) 
1 T(k s , -k s + k t ,k' -k' - kA , 

9 1 r T(k s ,k' s ,k,,-k s -k' s -k ; ) 1 

25 9NL = 18 ft, PMPWPfr) 3 T7VL (R21) 

where in the expression for g^L we fix |k s | = |k' s | = |k s + k' s \. Substituting the power spectrum, 
bispectrum and trispectrum into Eq. (B.19)-(B.21) gives 



If 

5 JNL 2(/(1)) 



7 (2) J / (4) /ggV, au g^ F ! / /( 2) /( 3) (/(2)) 3\ 

/W) 2 I 4(/W) 2 2(/(!)) 3 i 4(/W) 3 Cg1 " 4 1 (/(i))3 (/(i))4 J ^oK^n--- 



(B.22) 



A / (3) , fl / (5) l(/ (3) ) 2 \ ^ 2 v , f (/ (3) ) 2 (/ (2) ) 2 / (3) , / (2) / (4) ^ T , 

25 9NL 6(/W) 3 ^12 (/(1))3 4 (/(i))4 J + ^6(/«) 4 S^ 1 )) 5 6(/(!)) 4 J Cg1 s) 

( 1 / (2) / (4) (/ (3) ) 2 , 1 (/ (2) ) 2 / (3) \ , , M , i Cggg. ^x , m «,x 

+ ll2(/(l))4 6(/W) 4 12 (/(1))5 (/(1))5 + {U.A6) 
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and 



r - (/(2))2 4- ( /(2)/(4) 2 (/(2))2/(3 M f 1 ( /(3) ) 2 + Iff r (k WB 24^ 

T7Vi - (7^7 + ^77«F" " (/d))5 J (Cg> + ^2 P7 + pjs J /cg( ^ (b - 24) 

/ (/(2))2 / (3) _ (/ (2))4 \ J Co (fe.)+J CG (^) 

V (/ (1) ) 5 (/ (1) ) 6 / 2 

The observed values of the non-Gaussian parameters /jvx, ^tvl, tjv£ include the scale-independent 
loop contributions at Q((C G )), rewriting the series Cjvg(x) = E„ ^/ (n) CG( x ) as Ctvg(x) = E n hM H u(Cg(^)), 
where H n are the Hermite polynomials (as we have done in Eq. (3.1)) cancels these lowest order 
loop terms. 

If instead we have = and / (2) ^ O.For /W ^ then, 

^ L = feS)7^ + --- ™ = (/< 2 ')^(fc)/ Co(k ;) + --- <wi = o + ... < R25 > 

C Mapping Between Statistics in Vl and Statistics in Vs 

We would like to calculate the local statistics, that is correlation functions of CtvgIs- First, we 
rewrite the locally observed non-Gaussian field as 

GVG(X)| 5 = (/« + f^CG,L + ^ (C 2 G ,L ~ (C G ,l)) + ^CkL + ■ ■ ■ J CG,S (C.l) 
+ ^ (/ (2) + + ^ " (Cg,l)) + ^C G ,L + • • • ) (C G ,5 " <C G ,S» 

+ ^ (V 3) + / (4) Cgl + ^ (Cg,l - <Cg,z» + • • • ) C C 3 

+ I [f {4) + f {5) <G,L + ^ (C G ,L " (CG,L)) + • • • ) (C G ,S " 3(C G ,5) 2 ) 
+ ... 

n( 2 ) »(3) n( 4 ) 

= S (1) Cg,s + /y 2 (Cg,s - <C G ,s» + + V (A* ~ + ■ • • (°- 2 ) 

where 5 = /'(Cg,l)- The coefficients g^- 1 are equal to /^"^ up to corrections 0(/^ n+1 ^CG,L)- So, the 
local statistics are similar to the global ones as long as the amplitude of the background mode obeys 
C^ L < m \f( n )/f( n + m ). 

Under the approximation that the coefficients are constant across the volume Vs we can 
use the expressions from §B with f g and I( G (k) — > s (k) : 

3 f 1 - 9i2) 4. ( ^ V 2 v 4. 9(2)5(3) / fc ^ 

1 ( g (2))3\ 
+ 4 I ( 5 (D)3 (5 (D)4 J^. S ^) + ... 



-3 

->G,S 
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1 ff(3) /l g(5) l (ff( 3 )) 2 \ 2 / (g (3))2 ( g (2))2 g (3) g (2) g (4) \ 

25 5JVL ' 5 6( 9 W) 3 + ^12 ( 5 (i))3 4 J (Cg > 5) + y6(gW)* 3(gW) 5 + e^ 1 )) 4 ) 

/ 1 gfflgffl (g (3))2 ! (g (2) ) 2 g (3) \ 1 (g(2))2 g (3) 

+ l v 12( 5 (i))4 6( 5 W) 4 + 12 ( 5 (i))5 J J fe.5^J + 6 (5 (i ))5 PCg,sW + ~- 

and 

/ (gfflggg (g(2))4 \ I CGs (k s ) + I CGS (k' s ) 
1/3 

where fc s , fc/ > 2tt/V s are the long and short wavelength modes used to measure /nl, 9nl tjvl 
within V5. 
So that 

+ ylfi?) ~ 70T J ~ ^Gt,l)) + ^|(iy - ^(2y^) ~ ^,5)) 



and when /W = 0, / 0, 

fNL\s _ Z<G,L + 2I C G ,s(h) + 



fNL f {2) {kc,s+2C 2 G , L ){kc,s+2t 2 G,L 



D Diagrammatic Rules for Fourier-Space Expressions 

The real-space diagrams shown in Appendix A can also be calculated in Fourier space: 

f " d 3 ki 

real-space expression of (xi, . . . x n ) = J J^J ^ ^ [fe-space expression] e 1 ^^'^ . (D.l) 

Here we will show how momentum-space expressions such as those given above in Appendix B for 
the bispectrum and trispectrum can be quickly recovered from their corresponding diagrams. 

For the local ansatz £tvg( x ) = Ylm ^1/ (Cg ( x ) ~~ (Cg ( x )))' a particular n-point function is 
given by 

f(mi) f(m 2 ) f{m„) 

(Ck 1 c k2 ---co = EE---E mi , m , ■•• ^((CG x ) kl (CG 2 ) k2 . . . «s»)o, (d.2) 



mi 7712 mrj 



- 30 - 



where the modes are given by the convolution integrals 

„ m— 1 73 m—1 

(c") k = y n (^c G ( P .)c G (k- 5: Pi ). (d.3) 

A given term in (D.2) is specified (up to permutations in the kj) by a set of numbers V m , where V m 
is the number of times the term appears. Thus, Yl Vm = n - Restricting to the connected part 
(CkiCk 2 • • • Ck„)c imposes the condition \ ^ mV m — n + 1 > 0. This is not a sufficient condition for 
the contribution to be connected; the contractions must be made so that the corresponding diagram 
is connected, as discussed below. Note that ^ mV m is even for nonzero contributions. 

There are m — 1 momentum-space convolution integrals for each of the (Cg )k^ f° r a total of 
^2mV m — n integrals. The factors of Qg m (D.2) are contracted using Wick's theorem, giving 
\ S m Vm — 1 delta functions, not counting the final overall <5 3 (^] kj), so in the final expression there 
are L = ^ ^ rnV m — n + 1 integrals remaining. This is the number of loops that will appear in the 
diagram. If L = the graph is a tree graph. Tree graphs dominate contributions to the n-point 
functions for a weakly non-Gaussian series. If L < the diagram is disconnected. 

The (Cg )ki factors will be represented as m-point vertices, with V m of each type in the diagram, 
and a total of n vertices. V\ = E denotes the number times the linear term contributes; these 1- 
point vertices appear as external lines. Finally, each contraction between two factors of C,g yields a 
factor of the power spectrum and is represented by a line connecting two vertices, with a total of 
P = 2 m Vm = n + L — 1 lines. 

The rules for diagrams are as follows: 

1. Assign a momentum label kj (i = l,2...n) to each vertex, including external 1-point vertices. 
Each m-point vertex is equivalent to a factor of /( m ). (The ^ is cancelled by the m! ways of 
contracting into the vertex.) 

2. Assign a momentum label to each line, with a direction. Lines contracted with 1-point vertices 
share their momentum label. L internal lines can be labelled with integrated momenta pj 
(j = 1,2. ..£»); these can be chosen arbitrarily among the lines forming loops. The remaining 
P — E — L = (n — E) — 1 internal lines can be labelled with momenta k/ + ^ , where k/ 
is the momentum of one of the vertices contracting with the line (either can be chosen), and 
q/% {k = 1,2.. .m — 1) denote the incoming momenta of the other lines being contracted into 
that vertex. This imposes momentum conservation at each vertex. These labels can be made 
by working into the diagram starting from the external lines. Each line is then equivalent to 
a factor of Pc(q), where q is the momentum for that line. 

3. Integrate over the loop momenta by adding a factor f for each loop. Note that loops at 
a single vertex contribute a factor (Cg)- 

4. Divide by the symmetry factor of the diagram. As in standard quantum field theory, the 
symmetry factor is determined by counting the number of ways of exchanging identical vertices 
or identical lines, as well as lines contracted at a single vertex. 
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5. Sum over permutations of the k, (momenta for the vertices). Sum over connected diagrams 
with n vertices, to desired loop order, or level of approximation, and multiply by (27r) 3 e) 3 (^ kj) 
to obtain the n-point function (CkiCk 2 • ■ • Ck n )c- 

These diagrams are essentially equivalent to those considered in [49], where the more general 
case of multiple fields contributing to the curvature perturbation was considered. 
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